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The n-point function for the integral over unitary matrices with Itzykson-Zuber mea- 
sure is reduced to the integral over Gelfand-Tzetlin table; integrand (for generic n) is given 
by linear exponential times rational function. For n=2 and in some cases for n > 2 later 
in fact is the polynomial and this allows to give an explicit and simple expression for all 
2-point and a set of n-point functions. For the most general n-point function a simple 
linear differential equation is constructed. 
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1. Introduction 



In this letter I'll consider the following correlation function : 



< 9^ljl9th■■■9^^j^9k„l^ >= J Wig)] exp[Tr{gMg+N)]g,^j^g+^^...g,^j^g+^i^. (1.1) 

Here g is the N dimensional unitary matrix and M and N are Hermitian. Measure of 
integration is Haar measure. Without lack of generality we could assume that M and N 
are diagonal. 

For the case of n=0 (partition function) this integral was calculated by Harish-Chandra 
1^ and Itzykson and Zuber [|2| long time ago. Here we will use the method previously 
used in similar problem in ; this simple algebraic method is known in literature on 
representation theory since 1950 ||^]. Let me mention that the main motivation to look on 
the integral (|1 . 1| ) is related to investigation of Kazakov-Migdal model ; also , this kind 
of integrals might be interesting for string theory related matrix models f^. 

I'll show that the integral ( p. . 1| ) can be written in the form: 

^ „N-lN-k N N-k+l N-k 

^ k=l i=l k=l i=l i=l 

(1.2) 

where A(M) is the Vandermond determinant constructed from eigenvalues of matrix M, ml, 
are Gelfand-Tzetlin coordinates defining the convex body (see fig.l) 

<>™^'>mfc+i> (1-3) 

with m° being the eigenvalues of matrix M (we assume that they are ordered: > m2 > ... > 
), Nk are eigenvalues of matrix A^, C is the number (see (2.12) ) and R[m] is given (in 
general case ) by rational function on GT table. This function will be described explicitly 
in section 3 together with linear differential equation for (|1.2| ) (I'll give the corresponding 
algorithm) . 

R[m] reduces to simple polynomial in the case of ii = ^2 = ... = ■in = 1 and this allows 
to give the explicit formula (we denote by A(A''2, ...,N ~ N) the Vandermond determinant 
for eigenvalues A^2, -^3, Nn): 
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+ + + _ r r r 

exp[(^m^ - X] m^)iVi]det[exp(m-iVj)] 

fc=l A:=l 

Here z = 1, 2, — 1; j = 2, 3, A^. All nonzero 2-point functions < gik9ki > a set 
of multi-point functions < 9iji9k„i---9ij„9k„i > are obtained from (|1.4| ) by permutation of 
eigenvalues after integration over in 



2. Gelfand-Tzetlin parametrization 

Let us denote by X the combination that enters in IZ measure: 

X = gMg+. 

X is the element of coadjoint orbit of unitary group; this orbit is labeled by eigenvalues of 
M : rrii > ... > . We will now introduce the coordinates on group and on G/H ,with H 
being the Cartan subgroup defined by M , so called Gelfand-Tzetlin coordinates (for the 
reason why we call this coordinates GTC see [§]). 

Let tti be the basis in N dimensional complex plane, , and be fixed basis , in 
which matrix M is diagonal with ordered eigenvalues m°. Than we can write 

N N 

M = ^m°eie+,X = ^m°aia+,^ji = {ei,aj). (2.1) 

i=l i=l 

o-i = fi'Ci, ( , ) is the scalar product: (ax,?/) = a*{x,y); {x,ay) = a{x,y). Haar measure 
could be parametrized by the vectors a^, dfj,{g) = dfi{ai, ...,aN)- Then 

(A^ - 1)! 

dfj,{ai, ...,aN) = —i—y^dti...dtN-id9i...d6Nd^{a2,...,aN), 

^ Integral in ( |l.4| ) is easy to calculate, but I don't think that integrated version looks simpler 
than this. 
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tk = \9ik?,k=l,...,N-l- (2.2) 



0k = aTggik,k= 1,....,N. 

One can move to the new variables mj from ti in following way: suppose is the 
eigenvalue of matrix PMP, where P is the projector on the subspace, orthogonal to vector 
ai. This means that if / is corresponding eigenvector, / = Ylk=i ^feefe,one has 

N 

PMf = mf, Mf = m^f + bai,f = J2 hek- 

k=l 

Taking the scalar product of Mf with //, and e/. we get 



N 



^9lk ,2 \9lk 



^fc = — 1'^ = 

my — m'- < ^ 



|2 



m\-m}-' ^ (mfc - vaf ' 

From the orthogonality condition, (/, ai) = 0, one obtains equation, which relates g\k and 
eigenvalues by 

I 1 2 

— 1 = 0)*fc = ^I/C 5 

AT 

5]tfc = l. (2.3) 

fe=i 

This equation describes one-to-one correspondence of coordinate system and 
m^; from it immediately follows that eigenvalues obey inequalities 

m°>mj>m°+i, (2.4) 
and this intervals are filled densely. Moreover after some simple algebra one finds relation 

= (2.5) 

For to move to coordinates m}- we also need the Jacobian for integration measure. 
From above expression it could be easily obtained and finally we have 
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(2.6) 



(27r)^ 

In addition one can show also that foUowing formula holds 



(i/z(ai, Cat) = — - — -jj—J{t,m^)dm\...dm]^_id9i...d6]\[dn{a2,---,a]y). 



N N N N-1 

^11 = ^gikrrikgt^ = ^tkrril = - ^ (2.7) 

i=l 1=1 i=l 1=1 

This completes the first step in descent procedure, which embeds the A^— 1 dimensional 
unitary group in N dimensional one. After this step we have — 1 dimensional space 
spanned by ortho normal vectors a2, oat, diagonal matrix = PMP with eigenvalues 
m\, rn\i_i and eigenvectors / . It is simple to continue the descent further; let me give 
the necessary expressions after the descent procedure is completed down to the level i 



N-l N-2 N-l+l 

'32' 



1 1 1 

(2. 



[jq] 9 = Jq Jq + 1 



^q^fk ^) = VtkqexpiiOqk), (2.9) 



( fq fg-l^ _ „,q (^9' fk ) 

\J i ■> J k ) "i q-\ q ' 



_k ; 

ml ' - ml 

where by a and t we have denoted following products 



]\^=i i.'m] - ml ^) 



^ one can derive all this expressions using the procedure described above, or (after some mod- 
ifications) they could be extracted from see also @ 



4 



The variables are (as it follows from descent procedure, see ( p.4|) ) from Gelfand-Tzetlin 
table 

m\ > m^+i > ml, (2.11) 

and 9 are just angle variables, 9f = [0, 2tt]. Using all this one could check by straightforward 
calculation that 



N N N-k+1 N-k 

1=1 1=1 i=l 1=1 

^ ^ N-lN-k N N-k+1 

= TT-idii^Aim n n ^< n n 

l^^J 2 ^ ^ k=l i=l k=l i=l 

Cn = (iV- 1)!(A^-2)!...2!1!. 

The meaning of symbols fl is simple; f.e. fl forms the orthonormal basis in the space 
spanned by eigenvectors of matrix = diag{rn}f^) 



3. Correlation Functions 



Having at the hand the expressions (|27^) , (|2 . 9| ) , (|2 . 1 0| ) , (|2 . 1 2| ) it is easy to prove the 



statements about ( |1 . 1|) , announced in introduction, i.e. derive explicit formulas for all 
2-point functions and for a set with n > 2 ( |1.4| ) and present the integral representation 
of general n-point function in terms of GT table (|1.2|) . Also one could give the linear 
differential equation for general n-point function . But first I will start with 0-point function 
and will give the simple derivation of IZ formula: i 



< 1 >= y" [(i//(^)]exp[tr(^M^+iV)]; (3.I) 



integration is over Gelfand-Tzetlin table (|1.3| ). The result follows immediately from the 
facts we have learned in previous section (we need just (|2.12|) ) and from two simple obser- 
vations : i. Integrand doesn't depends on angle variables 6, so it is the linear exponential 
on GT table: 



^ this is one more derivation of well-known IZ formula but I decide to include it here for 
completeness of presentation ; also this kind of derivation might be itself useful, i.e. for large N 
limit. 
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< 1 >= 



„ ^ N N-k+1 N-k 



A(M) 

ii.the integration over the bottom coordinate m^of 



J n [dm] exp[X;( - E ^')^^]- (3.2) 



GT k=l i=l i=l 



exp([(mf-^ + m^-^) - m^] A^jv-i)exp(mjviVjv) 



leads to 

det [exp (mf^ ~ A^iv_ j + 1 ^ 



Nn - Nn-1 

and variables m are from the second (counting from bottom) line;z, j = 1, 2. Now, according 
(^.2|), we have to multiply this expression on exponential 



3 2 



exp(5^mf-^-5^mf-i)7V^_2 



1=1 1=1 

and integrate over rn^~^ ; once again we get the similar result: 

det[exp(mf-^Arjv_,+i)] 

with i,j = 1,2,3. Thus, using the induction method one easily proves that following 
formula holds: 

'"^ N-z+ii _rV^_7V-i ^^N-i+ii.r det[exp(mf-'+^A^Ar_j+i)] _ 



1=1 i=l 

det [exp ( rrip ~ ^ A^at _ q_|_ 



A{Nn,..,Nn-i+i) 

(3.3) 

here p,? = 1, ...J^iJ = 1, (l - 1). 

This completes the derivation of IZ formula in our approach, based on GT parametriza- 

tion: 

where we have denoted by Mi. 

Let us now move to discussion of correlation functions. First about general properties: 
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All nonzero correlation functions in (|1.1| ) have jm = km and (/i, = P{ii, i^), where 
P is the element of permutation group; of course the "reversed" statement , when = Im 
and (fci, kn) = -P'Oi, ■■■,jn) , is equivalent to above. 

To demonstrate this simple fact let us multiply right hand side in ([Ll])by 



where Hi, H2 are from U{1)^ . Because the IZ measure is invariant under the transforma- 
tion g' = HigH2 we can perform this transformation in the integrand ; thus each operator 
inside the correlation function will be multiplied by Hi from left and H2 from right (her- 
mitian conjugates will be changed correspondingly). Because the integration over Hi, H2 
is over U{1) angles we see that the integral will be always zero, except the cases that stated 
above. 

The expressions derived at the end of section 2 and the statement ( |3 . 3| ) immediately 
leads to prove of (|1.4D . We have 



GT 

n'^-Vml - 77)0 1 ^ N-k+1 N-k 

■n.^,„K-0 ^tt^ ^ t,'^ A(M)A(iV„...,7V^) 
j^^ ^^1 YifJi (mj - TO?, ) U^=~i\rnl - m% ) ]\iSi{m] - mj^ ) 



[m^ — mP- ) 



' Y\i^,M - ml) nz/,.(™? - o '"Yii^jM^i - m] 



N N-1 



exp[(^m^ - ml)Ni]det[exp{mlNj)]5j^k,-Sj„k„ 
k=i k=i 

(3.5) 

Here i = 1,2, .., N — 1; , j = 2, N. As I have mentioned in introduction the last integral 
is easy to calculate; it is clear form integral representation that the result is the sum of 
quadratic exponentials of eigenvalues Mi , Ni times rational functions and it is regular when 
two eigenvalues of M coincide.!^ So we could permute the eigenvalues of M to obtain the all 
other correlation functions < \gik\'^ >, as well as all n-point functions with ii = 12 = ••• = in 



^ we have used the coordinates tied with matrix M but we could use in the same integral 
the GT coordinates related to matrix A'^; result shouldn't depend on this choice; note that GT 
coordinates for M eigenvalues of A'^ aren't ordered. 
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(and others that are related to later by interchanging matrix M with matrix A^) after last 



integration in (|3.5|) . 

Now we will discuss the most general n-point function, n > 2, the one that couldn't be 
reduced to (|3.5| ) , and will prove ( |1.2|) . Using general property discussed above we will order 
the correlation function ( |l.l|) by the set of indices (zi, > ^2 > ••• > 'in- According 

( ^78|) we will draw the two copies of Yung type diagram : one for set (ii, ji, jVi) 
another for (zi, i^, ki, kn) (see fig. 2). Each line in the diagram has length i^, and thus 
we write numbers (h, in) in the boxes at the right end of each line in each diagram. In 
the first diagram we put numbers (ji, jVi) in most left boxes (again from up to down, 
but now the numbers aren't necessarily ordered); equivalently in the second diagram we 
put at the same place numbers {ki, kn) = -P(ii, ■■■,jn)- So each line in first diagram has 
iq in the right and jq in the left, at the same time the same line on second diagram has 
iq at the right and P{jq) at the left. Now we will fill all other boxes in the first diagram 
according following procedure : pick up all lines with length in (some of the numbers from 
set i might be equal, so we will have several lines with the same length) and put next from 
in the numbers {p}^_i, ...,pf^"i) from up to down with requirement that 

1 <Pj„-1 <N-in + l. 

Ai^ is the number of lines of length z^.We continue this procedure before reaching the 
point with in — P = ^n-i, where p counts the boxes from right to left on the top line of 
diagram. After this the height of each column is larger, so next column we will fill with 

A ■ 

numbers (p^„_^_i, ...,PiJ"i^), with Ai^_-^ being the number of lines of length in plus the 
number of lines of length i^_i.This procedure we will continue to the end, when all boxes 
in Yung diagram are filled. The numbers p should satisfy unequality 

l<pf_,<N-q + l. 

Let us call this diagram Y'^^{p). The same procedure we will repeat with the second dia- 
gram, assigning the numbers p' to empty boxes and defining V^'^^^^p'). 
This diagrams are just convenient representation of products 

9llh-9tr.jr.^ (3.6) 



when they are viewed in terms of sum in ( p.8|) . For to proceed further we have to assign 
to each box in the diagram the factor 



T .+1 (3-8) 

Jq Jq + i 

and then multiply over q in each line; the same we should do with the second diagram, 
but assigning to the boxes the complex conjugate of (|3.8| ). So, for each diagram we have 
assigned the particular product dictated by above procedure and the diagram. Let us 
call these factors F*-^(y) and F'^^^^\Y'). Finally, for to obtain the factor R coming from 
insertion of (|3.6| ) and (|3.7|) in IZ integral we have to sum over all possible diagrams the 
objects F{Y) and F{Y') ((with fixed (zi, zat); (ji, jat); P ) and then multiply: 

i? = 5^F(y)F(r)- (3.9) 

YY' 

The described procedure looks very complicated but result of integration over angle 
variables makes it simpler. Let us take the product of Fi and F2 with some particular 
set of p and p' and integrate over angle variables with measure ( |2.5| ). Because we already 
know that the dependance on angles is coming only from pre-exponent in ( p.. 1|) (see ( |2.5| ) 
and ( |3.2| )), this is the only angle integral we have to calculate. From other side the angle 
dependance in (|3.8[) is due to (|2.10[) , thus we have simple angle integral to calculate: fix 



the column (let say 1-th) in both diagrams and denote by G sum over this column of the 
angles ^p;,z; denote by 0' same sum for second diagram. The angle integral that enters in 
(Ol) ^^^^ column is just 

n (fi%xp(e-e'), 

'•ow [Z] 

which is zero if the set p' on the second diagram in given column (1-th in this case) is not 
related to set of p in this first one by permutation 

(p^...,p;^0 = p(p,\...,pfO- 

Later relates the numbers in second diagram to numbers in first by simple law: we 
should consider only those pair of diagrams, when numbers in fixed column of second are 
related to numbers in the same column of first diagram by some element of permutation 
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group, so for each diagram Y we have set of mappings r defined by above law: Y' = t{Y). 
Thus we have for R after integration over aU angle variables: 



Rq = Ro[m] =J2F^'''\y]Fo'''^'^\r{Y)]. (3.10) 

Here subscript means that we have to put to zero all angles in ( |3.8| ) ; after this R becomes 
the function on GT table, Rq [m] Now it is easy to see that Rq [m] is the rational function. 
Simply after applying the mapping r to y there are only moduli squares of (^]^) and ( |2.1Q ) 
left in( |3.10| ); this means that -Ro[to] is rational. On this we complete the construction of 
R. Going back to (|l.l|) in parametrization ( |2.12| ) we have (for nonzero correlators): 



Cn f ^ 
< 9^lJl9tpiJ,y■■9^„Jr.9tp(^J„) >= J l[dmRo[m]exp[J2XkkNki 



C 



N 



A(M) 



GT k=l 
„N-lN-k N-1 N-k+1 N-k 

fc=l i=l A;=0 1=1 i=l 

(3.11) 

and here indices for R[m] are arranged according for those in right hand side. Thus we 
have derived the formula announced in the introduction. 

The integrand in (|3.11| ) still looks complicated, (even having the explicit algorithm 
of its construction), except the one when ii = i2 = ■■ = in', but it seems that it could be 
simplified using algebraic identities for particular rational functions that enter in Rq I. It 
would be also interesting to have some geometric interpretation of above construction. At 
the moment best we could extract from ( p.ll|) in general case is following: first we could 
extract from i?o[™] two polynomials Pi and P2 by defining 

r T Pi [ml 

«oM = ^ (3.12) 

( ^.12| ) uniquely defines Pi, P2; there is a simple formula for P2 in the case of general n-point 
function P2 

P2 = n Il« - mf^T" ll« - ^1)'" lliK'^ - mr^Y'^ (3.13) 

q] i^j k^l p^t 



^ this should be viewed as conjecture based on knowledge n-point function ( |3.5| ); i.e in our basis 
the expression for < \gik\'^ > has similar structure as general correlation function in ( p.ll| ) and at 
the same time answer obtained from the one with z = 1 by permutation of eigenvalues of M is 
simple. 



10 



where a set [q] and numbers aq,bq,Cq are fixed by the particular correlator under consid- 
eration. Then consider more general then (|3.11|) inegral: introduce for each point on GT 
table corresponding "dual" coordinate and define / by 



7V-1 N-k 



I{z) = I \[dmR[m]e^^[Y^ z^m\] (3.14) 

GT k=l i=l 

From I{z) integral ( |3.11 ) is obtained by requirement that all coordinate on the "z" table 
in given row k are equal to A^^ — Nk-i; From other side I obeys the linear differential 
equation 

P2{d,)I = Pi{d,)Io, lo = I{R = 1), (3.15) 

so that right hand side in ( |3.15| )is known function. Thus, some particular solution of 
( p.l5|) for the values of z , that I have defined previously, should coincide with ( p.ll|) and 
therefore with ( |1 . IJ ) . 



4. Conclusion. 



One should note that the expressions (|1.4[ ),( [L2D are very simple. This simplicity of 



course should be related to the fact that the integral ( |1.1|) actually is not over unitary 
group G, but rather over it's left/right coset. We have used only part of symmetries when 
we tied the coordinate system to matrix M (and not to matrix A^). Using the fact that 
result shouldn't depend on this choice additional identities could be obtained. It might be 
possible that later allows to get all nontrivial n-point functions for n > 2 from those that 
we have calculated (see ( |1.4| )). 

The simplicity of (|1.4|) , (|1 . 2| ) might be helpful for the applications, mainly for large 
limit; later is important both in string theory and Kazakov-Migdal model. In this 
respect I would like to note the elegancy of (|3.3|); it might lead to interesting equation in 
large limit. This results might be also helpful in the statistical models , discussed in0. 
A.Migdal have derived explicit formula for general 2-point function in the leading order 
of large A^ approximation in his approach based on Riemann-Hilbert method Q.I know 
from A.Morozov that he have conjectured simple formula for 2-point function that looks 
similar to the one obtained from the (|1.4| ) after integration over in the case of 2-point 
function. It would be interesting to find the geometric approach that could explain in more 
invariant way the results obtained in IZ model both for finite and large A^. 
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